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1 Introduction 



Understanding flow patterns in a convergent-divergent nozzle (the so called 
de Laval nozzle in engineering) is a prominent issue in aerodynamics and par- 
tial differential equations due to their numerous applications in practice, and 
closely connection with many difficult mathematical problems, such as mixed 
type equations and free boundary problems Pll3lll4|17 j. Since these ffow pat- 
terns are genuinely nonlinear, various physically significant special solutions 
of the corresponding mathematical problems play an important role in the 
theoretical analysis. For example, for the nearly spherical symmetric tran- 
sonic shocks and transonic shocks in a slowly varying nozzle, there are works 
of Chen et.al. [T1l5|15ll22p^ on potential fiow equation and complete Euler 
system based upon in essence two classes of special solutions [25J. There are 
also many progresses in the analysis of subsonic nozzle fiows, see, for instance, 
p!6|l2T] and reference therein. However, since no simple and physical special 
transonic-fiow solution is available, presently the study of subsonic-supersonic 
transonic fiow mainly utilized the methods of compensated compactness (see 
[51IHI17II21] and references therein). 

In [23], motivated by a significant work of Sibners [20], Yuan constructed var- 
ious interesting special solutions in a two-dimensional Riemannian manifold 
with "convergent-divergent" metric, which may be regarded as an approxi- 
mation of the de Laval nozzle. This paper also studied several boundary value 
problems of subsonic fiows in such a manifold. In the present paper we will 
further investigate the subsonic-supersonic fiow via the potential fiow equa- 
tion. We show that the special subsonic-supersonic transonic fiows are stable 
with respect to the small perturbation of the velocity potential function (i.e., 
the tangential velocity) at the entry (see Theorem [3]). 

The potential flow equation is a second order equation of elliptic-hyperbohc 
mixed type for transonic flow. For such equations, presently one of the main 
tool is the theory of positive symmetric systems and techniques of energy 
estimates, see, for instance, [8]|13fl4j and references therein. In this paper we 
employ the theory developed in [13] to show the solvability of linear problem, 
and then a nonlinear iteration argument to solve the nonlinear problem. 

We remark that in recent years there are many breakthroughs on partial dif- 
ferential equations of mixed type and degenerate elliptic type arising in dif- 
ferential geometry and physics, see, for example, [2!5) MTU|11II12II14|17] . For 
earlier developments in this fleld, one may also consult [T7|18lll9j and refer- 
ences therein. 

The rest of the paper is organized as follows. In Section 2 we formulate the 
problem, and study the properties of the coefficients of the potential flow 
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equations in the manifold. In Section 3 we solve the linear problem, and finally 
in Section 4 we state the main result, Theorem [5] and prove it. 



2 Formulation of the problem 



Let be the standard unit circle in R^, and Ai be the Riemannian manifold 
{{x^^x'^) G [—1,1] X S^} with a metric G = Qijdx^ ® dx^ = dx^ ® dx^ + 
n{x^Ydx'^ ® dx"^ . Here n{t) is a positive smooth function on [—1, 1] satisfies: 
(1) n"{t) > 0; (2) n'{t) < for t G (-1,0), n'{t) > for t G (0,1). Such 
a manifold Ai may be regarded as an approximation of a two-dimensional 
convergent-divergent nozzle, with Ai^ = fl {x^ ^ 0} respectively the 
divergent and convergent part. We also call S'^ = {k} x S^, k = —1,0,1 
respectively the entrance, throat and exit of Ai. Obviously dAi = U S^. 

Let p, p be functions in At represent respectively the pressure and density 
of gas fiow in Ai, and f be a vector field in At represent the velocity of 
the fiow. We consider polytropic and isentropic gas fiows, then p = Kp"' with 
K > 0, 7 > 1 two constants, and the speed of sound is c = ^/nqp^. Let v be 
the 1-form corresponding to v under the metric G. The fiow is irrotational if v 
is exact; That is, there exists a function i^'m. AA such that v = dip. Substituting 
this in the equation of conservation of mass div(pf ) = —d*{pv) = 0, where div 
and d* are respectively the divergence operator and codifferential operator in 
A4, then by the formula d*{pv) = pd*v — {dp, v), with (-, ■) the inner product 
of forms in A^, we have 



pAip = {dp,d<^), (1) 

where A = dd* + d*d is the Hodge Laplacian of forms. (Note that d*(f = 0.) 
By the Bernoulli's law which represents conservation of energy: 



{d^,dv) + ^p'^-' = Co, (2) 

where Cq is a positive constant, p in ([1]) can be expressed in terms of dip. So 
we may write ([T]) as a second order equation of ip. 

Indeed, in the {x^,x^) coordinates, we have 



Aip = -—d,{^g'^d,ip) 

= Tn \di{n{x^)diip) + -^^822^] ■ 

n{x^) \ n[x^) J 
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Here ^ = ^J det (gij), {g^^) is the inverse of ((/jj), and di = d^i, dij = didj. By 
differentiating ([2]) we have 



/I 

—dp = -[-diipdjipdkg'^ + dik^pdjcpg'^jdx^. 
p z 

Then by a straightforward calculation we obtain 



+n{x')n'{x') \^c' + -l^{d2y^fj d,^ = 0. 



(3) 



Direct computation yields that this equation is of elliptic type if the flow is 
subsonic (c^ > {diipY + {d2'^Y / n{x^Y) , and is of hyperbolic type if the flow 
is supersonic (c^ < (c^iv?)^ + {d2fY /n{x^Y). 

If the flow depends only on x^, then ([3]) indicates that iph = (pbi^^) satisfies 
the equation 

nix'){cl - idiipbf)du^b + n'ix')cldiipt = 0, (4) 

where Cb is the sonic speed corresponding to (pb- It can be shown that there 
are special flows (pb G C^{Ai) which are subsonic in and supersonic in 
A4~^, and dipb > 0, dupb > in (see [23j). We call such flows background 
solutions. The aim of this paper is to study stability of certain background 
solutions under perturbations of p on the entrance of Ai. 

Let p = p — Pb- Bj subtracting (jlj) from ([3]), we have 



k{Dp)dnp + b{Dp)di2P + d22P - a{x'^)dip = f{Dp), (5) 
where 



a{x^) 



njx'ydnVb 4 + clidmf + (7 - m^pb)' 
dip>b 4 
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(9) 



njx'fdn^, 4 + cljdmy + (7 - 

We will investigate the following problem: 



di(p. 



Eq. ^ in M with = + V^fc, (10) 
ip = g{x^)oYi T,^^ with ||5'||//5(si) small. (11) 

This is a Dirichlet problem of an elliptic-hyperbolic mixed equation. 

For background solution let r = {diipb/cbY be the square of Mach number. 
By Bernoulli's law we may compute 



di{k{Difb)) = -[(7 + ly - 2r + 2] < 

r — i 



:i2) 



in Ai. Here we also note the fact that dnipb = {diipb)n' / [n ■ (r — 1)), and 
n'/(r — 1) > in Ai, especially at the throat S*^ we have n'/(r — 1) = 
nn" /('J + 1) by L'Hospital's rule in calculus. We also may write 



W|l + r + (-,-l)r^l^^_ 

r — 1 

so there holds 



2a - ldi{k{Diph)) > 6i > (13) 
in Ai for a fixed number 6i and all positive number /. In addition, 



TlTl T 

2a + di{k{D^b)) = -[(7 - 3)r + 4] > 52 > (14) 

r — i 

in Ai. We remark that the constants 6i, 82 depend only on the specific back- 
ground solution ifi,. 
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3 Solvability of linear problem 



Let k, b, a, a, f be functions in M = [—1, 1] x with fiat metric. In tliis 
section we investigate tlie solvability of the following linear problem: 

Lu := kdiiu + hdi2U + ad22U — a{x^)diu = / in M, (15) 
M = on E-^ (16) 

We use to denote the Sobolev space W^''^{M), and the corresponding norm 
is written as ||-||^-. Note that C^{M) denotes the usual space of j times con- 
tinuous differential functions in M . 

Theorem 1 Let k > on S"^ and k < on S\ and k,b,a,a G C^{M). 
Suppose f G H'^{M), < s < 3 and there is a positive constant 6 such that 

a>5>0, dia < -6 <Oin M, (17) 
2a - (2j9 - l)dik > S > m M for p = 0,1, s. (18) 

Then there is a u > such that if 
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< u, \\d2a\\c, <u m M, (19) 



then there exists uniquely one solution u G to problem (|T5l) and f|T6|) and 
there holds the estimate 



ll«IUi<c||/IL- (20) 

To prove this, we follow the celebrated ideas presented in [T3p^ . That is, one 
first solves a boundary value problem of a mixed type equation which is elliptic 
at both the entry and exit. Then the above theorem can be demonstrated by 
an extension technique. 

Theorem 2 Under the assumptions of Theorem [H, but supposing k > on 
both and T,^ , then the following problem 



Lu = f in M, (21) 

u = on S-\ (22) 

diu = 0on (23) 

has uniquely one solution u G and it also satisfies (I2U]) . 



6 



Proof of TheoremUl Step 1. Uniqueness. Let / = 0. Multiplying to (fT5l) by 
diu and integrating the expression in [—1, 1] x S^, Note that = on 
we have 



/ 


(a + 


Jm 








-L 




<o 





dx^ — 



dx 



(24) 



via the integration by parts and divergence theorem. Since dia < —5 < in 
M., we infer that, by choosing v = 6/4, 

{2a + dik + d2b){diuf - dia{d2uf + 2d2adiud2U 
>{6- 2u){diu)' + {6- u){d2uy > ^-muf + {d2u)'). 

Hence we have 



(6/2) / [{d.uY + {d2uy] dx^dx' < 



M 



by (EID and therefore w = in M due to (|T6|) . 



(25) 



Step 2. Existence. Let M* = [-1,2] x We may extend k, 6, a, a, / to M* 
such that they still satisfy (fT8!) - (fT9l) and other assumptions in Theorem [21 
especially A; > on the new exit = {2} x S^, and 



H^{M*) 



< c 



Denoting the obtained extended operator in M* as L again, we consider the 
problem (1211) fl22]) together with boundary condition diU = on T?. By The- 



orem [21 there is a unique solution u* and < C 



. Obviously 



u 



u*\m is also a solution to problem f[T5|) and f[T6l) . This finishes the proof. 



The proof of Theorem [2] also follows in a similar way of [H] (Theorem 1.1, 
pp. 9-18), but needs some modifications to deal with the mixed derivative term 
bdi2U, the non-divergence term ad22U, and no lateral boundary in our case. 
For completeness and convenience of the readers, we sketch out the proofs. 
Some of the details in the proof are important in the analysis of the nonlinear 
problems. 

Proof of Theorem Step 1. Uniqueness. This may be proved by a similar 
method as in deriving (p5l) . 

Step 2. Approximate problem. To show existence of a weak solution to 
problem (pT]) - (!23|) . as in [H], we consider the following singular perturbation 
problem: 
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L,u := Lu + eduiu = /, 1 > e > in M, (26) 
u{-l, x^) = diu{~l, x'^) = 0, (27) 
diu{l,x'^) = 0, (28) 

where G and — > / in as e — > 0. 

Step 2.1. Galerkin's method. To show existence of a solution to problem 
(!26l) - (!28!l . we employ the Galerkin's method of finite dimensional approx- 
imation. Let {Fj(x^)} {i = 1,2, ■■■) be a complete system in if^(S^) and 
orthogonal in L^(S^). We may also assume that each Fj(x^) is smooth. Set 



N 

u'^'^x', x') = Y: ^(x^)r.(a;2), iV = 1, 2, ■ ■ ■ . 

i=l 

The functions Xj^'" {i = 1, ■ ■ ■ , N) are to be determined by a boundary value 
problem of a system of third order ODEs: 



/jL,M^'^-/,]r,(x2)dx2 = 0, j = l,2,---,iV, (29) 

'^(-1) = (Xf ^)'(-l) = 0, (30) 
(Xf ^)'(l) = 0. (31) 



Equation fl29p can be written explicitly as 



N r 



6( /^^ dx')ixf'r' + j: [(/^^ kY^Y, dx') (X/ 

+ [ljbY!Y, - aY^Y^) rfx^) (Xf '^)' 

+ aYl'Y, dx^^ (Xf •^) = f,Yj{x^) dx\ (32) 

S'tej* ^.i.i. Uniqueness. Now we show the solution to problem fl29l) - fl3T|) is 
unique. Indeed, multiplying to ( l29l) by (X^^''^)', summing up for j from 1 to X 
and integrating with respect to x^ on [—1, 1], we have 



-/ [L,u^''diu^^']dx^dx'^ = - [ f,diu^''dx^dx^. 

JM Jm 

Writing m^'^ simply as w, then integrating by parts and using diw{—l,x^) = 
diw{l,x^) = 0, we obtain that 
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-2 / [L.wdiw] dx^dx^ = f [{2a + d^k + d2b){diwf - {dia){d2wf 

JM JM 

+2d2adiud2u] dx^dx^ + 2e / {duwf dx^dx^ + / \a{d2uf] dx^. 

JM JT.-^ 

Therefore by Holder inequality and Young inequality, due to (ITQll . we have 

e\\d,M\l+\\Dw\\l<Cml<C'\\f\\l. (33) 
Since w(— = 0, it follows that 



w{x^ ,x^) = w^x^ ,x^) — w^—l^x^) = J diw{t,x'^)dt. 
By Holder inequality, 



\wix\x^)\^ < j ^ l^dt j ^ \d,wit,x^)\^dt < 2 y ^ |9im;(x\x2)|W 
Hence, we have 



I'^^llo^ / \w{x^ , x'^)\'^ dx^dx'^ 

JM 



<4 / \diw{x\x^)\'^dx^dx^ = 4:\\diw\\l. (34) 

JM 

Now it easily follows from (!33|) and (IMI) that 



e\\duw\\l+\\w\\l<C\\f\\l (35) 

for a constant C independent of and e. This proves uniqueness of solution 
to problem (!29l) - (!3T]) . Note the above estimate also holds for the solution -u^ 
of problem (|26l)-(|28l). 

Step 2.1.2. Existence and regularity. Now by Fredholm alternative of boundary 
value problems of ODEs, we can infer that problem (|29l) - (j3Tl) has a solution 
M^'*" which satisfies ( |35ll . Indeed, we may write (l29l) -(l3T]) equivalently as a 
boundary value problem of a first order system of ODEs with 3n unknowns 



—^=AX + J^, (36) 
BX = 0. (37) 

Here g R^" are column vectors, and A,B are 3n x 3n matrices. Let 

'^1, - ■ ■ , '^3n be a set of 3n linearly independent solutions of the homogeneous 
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equation ^ = AX, and Xq be a special solution of the nonhomogeneous 
equation ^ = AX + JF, then to show existence, we need to find 3n numbers 
Ci, ■ ■ ■ , C3„ such that 

3n 

satisfies SA* = 0, or equivalently, 

{BX^,---, BX,rr){c,, C3„)* = -EXq. (38) 

This is a linear algebraic system, and it is well known that it is uniquely 
solvable if and only if for Xq = 0, then ci = ■ ■ ■ = C3„ = 0. However, this is 
guaranteed by uniqueness we proved in step 2.1.1. 

Note that all the coefficients in the left side of fl32p belong to and the right 
side of (El belongs to C\ so u^'' E C\ 

Step 2.2. Solution of problem (1261) - (!28|) . By the uniform estimate (1351) . there 
exists a subsequence {u^^'"} converges weakly in to a m*" G i?^ as j — oo, 
and diiu^^'^ converges weakly in to duu'^ G L^. We show that -u^ is a weak 
solution of problem fl2^-(l2g]). 

For Xj{^^) ^ C'^([— 1, 1]), by multiplying it to fl29|) . summing for j from 1 to 
N , and integrating with respect to on [—1,1], one has 

/ {[9i(A;x^)+52(6x^)+«x^-69nX^]9i«^'^+92(ax^)92«^'^+/.x'^}rfa;Mx2 = 

JM 

after integration by parts, where = Y^f=i XjYj. Supposing X strongly 

in H^, let ^ oo in the above equality, we have 

JM 

= (39) 

for all X ^ ^ H^. Therefore u*^ is a weak solution to fl26|) - fl28|) . 

Next we show satisfies ([57]) and Indeed, since H^{M) ^ L'^{dM), 

M'^''^(— 1, x^) = and m^'^ ^ u*^ weakly in indicate that u'^{—l,x'^) = 
0. Since e is fixed presently, and ^^^'^([-l, 1]; ^2(8^)) ^ C{[-1,1]; L\S^)) 
(see Theorem 2, p.286 of 0), ([35]) implies that ^im^'^ G C([-1, 1]; ^2(8^)). 
Therefore by ([30]) and ([3l]), we obtain that ([27]) and ([28]) hold. 

5*^61? 5. Existence of weak solution of problem (IT^ anc? ([TB]) . 
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Step 3.1. Weak solution. Now for e G (0, 1] we have constructed weak solutions 
of problem (!26|) - (!28l) which also satisfy (!35|) . Then there is a subsequence {m^^ } 
converges weakly in to a m G as tj — > 0. Obviously we have 



ll«lli<Cl|/llo. (40) 

We claim u is a weak solution of problem f[T^ and flTBl) . To verify this, let us 
take e = ^ in ([32]) for any x^H'^^Hl. Note that 



M 



<Ce\\f\L^Q, 



we get 



/ {[dx{kx) + d2{hx) + ax]diu + d2{ax)d2u + fx}dxMx'' = (41) 
Jm 

By approximation this holds for all x ^ -f^o- 

Step 3.2. Boundary conditions. The fact that u satisfies the boundary condi- 
tion M = on S"^ can be deduced by the same argument as in step 2.2. 

Next, let = M n {{x^,x'^) : 1 - o" < < 1}. Since A; > on by 
assumption, for rather small a, the operator L is elliptic in Mg^.. We claim 
that 

\\diU^\\HHMU<C\\f\\o. (42) 
If this is true, then clearly diu = on S^. 
Now we prove (gJl). Let G C°°([-l, 1]) such that 



. u I 0, -1 <x^ <1- 3a, 
, 1 — 2a < X < 1, 

where /i is a large positive constant such that di{arj) > on E^. 

Denote m^'^ by w. Multiplying (!29l) by 7]dii{Xf'''), summing up for j from 1 
to A^, and integrating the equality on [—1, 1], we have 
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' frjduw dx^dx'^ 

M 

{di -eri{diiwf + -di{r]a){d2w)'^ - -ari{diwf - ar]d2wdi2'w 



M 



+d2{ar]d2wduw) + {krj - ]^ediri){diiwY + br]duwdi2W + ^di{ar]){diwf 

—d2{ar])d2wdiiw — -dii{a'r]){d2wY + a^(f^i2'W^)^} dx^dx^ 
>[ ri[{k-i^-g){dnwy + {a-i^){di2wf]dx^dx'^-Crg\\w\\l. (44) 

JM 

Here we used the fact the diw = di2W = on S^, and to control the term 
— (l/2)e9i?7(9iiw)^. Hence by applying Young inequality to the left side of the 
above inequality and choosing u, g small enough, we have 



/ v[idnwf + id^2wy]dxV<C\\f\\l. (45) 

This implies by taking N ^ oo. 
Step 4- Regularity. 

Step Regularity in elliptic region. Let = M fl {(x^x^) : — 1 < X"^ < 
— 1 + cr}. Since > on we may choose o" > small such that the 

operator L is elliptic in M2a. By standard theory of elliptic equations (i.e., 
the interior estimate and the estimate near boundary of elliptic equations) we 
have 



\u\\uHM,.)<C{\\u\\,+ \\f\\,)<C\\f\\,. 
\'^\\h2{M'^^) - ^ ll/llo • 



(46) 
(47) 



Next we derive an estimate oi w = u^''^ similar to ()45jl near S ^. Let 
e(xi) gCo°°([-1,1]) satisfy 



0, x^ e [-1, -1 + ct/4] U [-1 + 2(T, 1], 

1, x^ e [-1 + a/2,-l + 3a/2]. 



(4J 



By multiplying ^dniXp') to ([29D, summing up for j from 1 to and then 
integrating on [—1, 1] with respect to x^, by the first equality in (14^ we have 
similarly 
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JM 



> (5/16) / ^[{dnwf + (duwf] dx^dx^ - C \\w 



JM 



1 ' 



hence 



Ijlidnwf + iduwf] dx'dx' < C WfWl 



(49) 



Step 4- 3- Regularity in mixed-type region. Let G C^([— 1, 1]) such that 

fO, -1 <x^ < -l + a/2, 
0<^9<1 and ^9(x^) = J 1, -1 + 3(t/2 < < 1 - 3(t/2, 

[o, l-a/2<x^<l. 

Multiplying ([29]) by Wui{Xp'), summing up for j from 1 to A^, and integrat- 
ing the equality on [—1, 1], we get 



Here we may estimate those terms involving -t}' , d" by the estimates (H5l) and 
(119!), since (supp C [-1 + a/2, -1 + 3or/2] U [1 - 3(t/2, 1 - a/2]. Therefore 



Hence by Holder inequality and Young inequality, the right side of (ISOll is 
bounded below by 



/ diif^)duwdx'dx^ = [ {{a-ldik)^{duwY -l-M'iduwY 




+di2{a'd)d2wdiiw — di{h'd)di2wdiiw + £^9(9111^)^ 




/ \^M\dnwf + -a^\di2wf] dx^dx^ < C\\f\\ 
JM 2 2 



(51) 



(5/4) / ^\D{^^w)\'dx'dx'-C{\\w\\l + \\f\\l) 

JM 



(52) 



Therefore, it is easy to get 
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with C independent of A^, e. Letting N ^ oo and e — ^ 0, we obtain that 

/ ^\D(diu)\'' dx^dx^ < 11/11? • 
By ( pi . fHB]) . flTTI) and the above inequahty, we obtain that 

\\diu\\,<C\\f\\,. 

Since a > 5 > in M, by we may estimate d22U. Hence we conclude that 
\\uh<C\\f\\,. 

Remark 3 We observe that in deriving the estimate, we used the assump- 
tion that D'^a, D'^a G L°° , but just required that Dk, Db G L°°. 

Step 5. Higher regularity. The regularity in (l20l) for s = 3, 4 in the elliptic 
region is obvious (also can be obtained by multiplier ^dfw,r]diW for integer 
A; = 4, 6), and in the mixed-type region can be obtained in the same fashion 
as those to (1451) by multiplier ddiW for integer k = 5,7. 

For example, for estimate, by the above energy estimate technique, we 
have 

2 ._1._2 



— / dii{f'd)diiiw dx^dx^ = / ed^diuwY dx^dx 

JM JM 



+ I {M - ldi{kd) + ]-d2{bd)]{dniwf - ldi{ad){du2wf} dx^dx" 

JM I I / 

+ / \d2{a'°&)d\\2wd\\\w — 2d\{b-&)d\Yiwd\\\w\dx^dx^ 

JM 

+ / (0^9)922 W^iiiW + 2^12 (a'i?)<9i2w9iiiW - 29ii(a^9)(9i2W<9ii2W 

JM 

\d\x{a-&)dxwdxx\W — dxv{k'&)dxxwdxx\W — dxv{^'&)dyiwd\\\w\ dx^dx^ 

\1 A^^Aji. 



9ii(a'i9)(9iiw) dx dx 

M 



> / ^^{^s_y_a)\{dxxxwY^{dxx2wr\}dx^dx^-CAff2- (53) 



M 



An important fact is that the terms with underlines involve only Db, so we 
may use ( |T9l) (and C^) to get the inequality. Therefore by choosing a 

small we get 



^{{pxxxwY + {dxx2wY\ dx^dx' < C 11/11^ . (54) 

M 

Combing an estimate obtained similarly in the elliptic domain, we have 
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\\duw\\,<C 



2 • 



(55) 



Letting ^ oo and e — 0, the weak limit u of the sequence {w = m 
which is a solution of (|2T|) . also satisfies 



ll^ll^illl < C 

By (EI]), we have 



(56) 



\\d22u\\,<C\\f\\, + C\\u\\, + C\\dnu\\,+ 
<C 



-duu 
a 









-0112U 


+ 


-di22U 


a 





a 



(57) 



for a constant C depending on 6. By (HM . we have ^di22U ^ < {i^/S) \\d22u\\^ < 
\\d22u\\, . Also, there holds -dii2U < C \\diiu\\.. So we get 

a Q 1 



l|522«lll<C„.„2 

and the regularity 



(5J 



l^^lls < ll^lla + + ||c^22^i||i < C 



Similar argument works for estimate. This finishes the proof. 



(59) 
□ 



Remark 4 We see here that we can not apply directly Theorem\^to ([5]) since 
dia < —6 < in f[T7l) does not hold in M. However, we need this assumption 
to control the mixed derivative term bdi2U. 

An important observation to TheoremUl is that its conditions are not invari- 
ant under multiplication of a positive function to the mixed type equation (fT5l) 
Therefore it is expected to find an appropriate multiplier to (j5]) such that The- 
orem [I] works. 

Now we choose h{x^) = e~^^^ with /i > a small constant (depending only 
on 62 and ||/i;||j;^oo), which is a bounded smooth positive function in A4, then 
it is easy to check that there is a positive constant 6* such that 



h'{x^) <-6*<0 in M, 

2h{x^)a{x^) + di{h{x^)k{D^b)) > 5* > Oin M, 
2ah - ldi{hk{Dipb)) > 5* > in M, 



(60) 
(61) 
(62) 
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where / > is less than a fixed number (say, / < 6). 

4 Solvability of nonlinear problem and main result 

Now we prove the following stability result of transonic potential fiows in Ai. 

Theorem 5 Let (ff, & he a background transonic flow in TVI . Then there 
exist positive constant C and Eq depending only on iph such that if 

ip = ipf, + g{x^) on (63) 

for any g G H^{S^) and \\g\\^ < £^ < £^0; then problem ([2]) and (!63|) has 
uniquely one solution with 

||<^-<^6||4<C^- (64) 

This result shows that the subsonic-supersonic fiow is stable under pertur- 
bations of the velocity potential function at the entry. More physically, since 
d2'p/n{—l) is the velocity of the fiow along direction, so gi^x"^) measures 
the perturbation of the tangential velocity along the entry. So we may claim 
that the special transonic fiow is stable under small variation of the tangential 
velocity at the entrance. 

PROOF. The proof is based on a nonlinear iteration scheme. 
Step 1. We define the iteration set as 

E, = {^e H\M) : 11^ - y^J^ <k< k^}, 

where kq is a small positive constant to be specified later. It is straightforward 
to check that there hold 

\\f{Dv)\\,<C,K\ (65) 



13 

/(D^(i)) - /(Z}^(2)) ^ < CoK ^^'^ - if^-'^ ^ (66) 



2 



(1) _ m(2) 



for any Lf, Lp^^\ ip^'^'> E E^. 
Step 2. Let 
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(h ■ M^){ip) := h{x')k{Dip)dii0 + h{x')b{Dip)di20 
+h{x^)d22'f — h{x^)a{x^)di(f 
= hix')f{Dcp), (67) 

where = (p — ipi,, and f{Dip) is defined by ffTOj) . By considering cj) = (p — 
g{x^) as the unknown, problem ([I])© and (JHSD is equivalent to the following 
problem: 



{hM^){^) = h{x^)f{D^) + {hM^){g{x^))m M, (68) 
= on (69) 

Here ip = (f) + g^x"^) + fb, = f — ft and note that 



\M^{g)\\^<Coe, 
M^w{g)-M^,,^{g) < CqE 



for 



(70) 
(71) 



Step 3. Now by ( l60|) -(l62l). we may choose kq so small such that for any ip e Ef,, 
there hold in M the following inequalities: 



h{x^) > 6*/2 > 0, (72) 

2h{x^)a{x^) + 9i(/i(x^)A;(D^)) > 572 > 0, (73) 

dih{x^) < -6*/2 < 0, (74) 

2ah - ldi{hk{Dip)) > 5* /2 > / = 0,l,---,5, (75) 

\\b{Dcp)\\^<i.* = 6*/A. (76) 

S'tej* 4- Then for any ip E Ef^, we solve the following linear problem of : 

(/iM^)(0) = h{x')f{D^) + {hM^){g{x^)) in M, (77) 

= on (78) 



By (j72l) -( |76l) and Theorem [H and the analysis in step 6 below, there exists 
uniquely one solution G and it satisfies 



Now choosing Sq < 1/(8Cq) and k = ACqE < kq (that is, C = 4Co), we get a 
(p = + g + ipb with ||(^||4 < K. Therefore we established a mapping T : ip ^ Lp 
on E^. 
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Step 5. Next we will show that T is contractive on in norm. 

Let e E^, T{ip^^) = (^W, and = ^^'^-g-^b, i = 1,2. Then 0^-0(2) 
satisfies the following problem 



hMn, (0(2))_/,M^,,, (0(2))] [/(D^(1))_/(DV9(^) 



= 



in M, 

on S"^ 



Note that 



M^(i)0(2) _M^(2) 



(79) 



then Theorem [T] and the results in step 6 below yield 



< Cis 



(80) 



By choosing £o further small, we obtain contraction. Then by a simple gener- 
alized Banach fixed point theorem, we proved Theorem [51 

Step 6. Solving H'^ solution of (177]) and (1751) wt/i /oti;er regular coefficients. 

For simplicity, we may write 

M^(0) := + Oi(D^)]9n0 + O2(/^<^)9i20 + ^220 - adi^ = F{Dip),{81) 

where k = k{Diph), Oi{D0) = k{Dip) - k{DLph), 02(00) = b{Dip), F = 
f{D(p) + M^g, and there holds 



Oillg <Ck< Ckq 



(82) 



for i = 1,2 and G -E^. Since h is smooth and bounded away from zero, the 
solvabihty of (IHI])(I7HD is equivalent to (177D(17gl). 

By Sobolev embedding theorem we have C^, so due to Remark [31 in 

Theorem m (jSTl) . (1751) has uniquely one solution and 



< C'llFll^ 



(83) 



However, we can not infer from Theorem [H directly the existence of if^ so- 
lutions and estimate like ^ < CUFUg. So we need the following apriori 
estimates and approximation argument. 
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Step 6.1. An apriori estimate. Now suppose there is a solution G to 
problem (IHTl) (ITHjl . we show that if kq is sufficiently small, then there holds the 
following estimate 



<C||F||^_^, forp = 3,4. 



The proof of the case p = 3 is similar to the following p = 4 case and even 
more simple (which should use the fact fl83l) ). Now we suppose fjH^ holds for 
p = 3 and to prove p = 4 case. 

First, for small a and kq, the operator is elliptic in M^cr- By the regularity 
theory of elliptic equation, we have 



< CIIFIL 



(85) 



Now let C G 1]) satisfy 



0<C<1 aud C(x') = r' 
~ [1, -l + 2a<x^<l. 

Denote u = ((j). Then it satisfies the following equation 



M^{u) = F:=CF + C'k(f) + C{2kdi(j) + bd2(j) - 
Since suppC' C (—1 + a, — 1 + 2a), by fl85l) . we have 



(86) 



F <C\\F\ 



3 • 



^7) 



Then by differentiating (15^ with respect to twice, we get w = duu satisfies 

M^p{w) + 2dikdiw = diiF + [—dukduu + duadiu + 2diadiiu\ 

-[^iiOi^iiM + dn02di2U + 2di0idiuu + 2di02du2u] 



as well as w = on S ^ by the cut-off function (. 

We note that the operator + 2dikdi in the left hand side of the above 
equation also satisfies the assumptions of Theorem[T](by multiplying a suitable 
positive function). So we have 
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Ikll2<c^ll^lli- 

We now estimate F term by term. 
(1) Obviously by (IHTj) we have 



9nF 



< 



<cm^ 



(2) Since /c, a G C^, we get 



-diikdiiu + diiadiu + 2diadiiu ^ < C \u\^ < C \\F\\^ . 
Here we use ( 183!) and (!84l) = 3) for the second inequahty. 
(3) We first recall the inequahty 



||ut>||-,^ < C \\u\\^ 
provided m G /f^ v e (see [14j, p.73). Then by (|82|) 

||5iiOi5nM||i<C||Oi||3||M||4<C/€||M||4, 

||ail02ai2w|li<C^||02||3 11^114 <C/€||m||4, 

||aiOianiM||i<C||Oi||3||n||4<CK||Mi|4, 
\\d,02dn2u\\,<C\\02\\^\\u\\,<CK\\u\\^. 

In all, we get ||F||i < C(||F||3 + k \\u\\^), hence 

||5n.i|l2<C^(l|i^ll3 + '^ll^ll4)- 
By ( I86I1 . we may then estimate 

||522«||2 < C(||F||3 + K \\d,^u\\^ + K \\U\\,) < Ci\\F\\3 + ^ 

So there holds 



II^IL < C{\\diiu\\^ + \\d22u\\) + ||^i||3 < C(||F||3 + K \\U\1). 

By choosing kq small, we can deduce that ||'u||4<C||F||3. Combing 
(l85i) . we can get (18^ for the case p = 4. 
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Step 6.2. Existence of solution by approximation. 

For (f G E^, Oi G i = 1,2, we approximate O, by {Of ^j^^ C C*^ such that 
of'* Oi{Dip) strongly in (so (|82|) holds). By Theorem [T5| the problem 



M(')(0) := [k + oP]9n0 + O^')9i20 + ^220 - ad,^ = F{D^), 
= on S"^ 

has uniquely one solution 0^'^ G H^. By the apriori estimate (18^ we have 
0(0 ^ < C* ll^llg for C independent of /. So there is a weak limit G 

of this sequence of approximate solutions. Then clearly is a H"^ solution of 
( 18T1) and fl78|) . and by Theorem [H this solution is unique. 

The proof of Theorem [5] is then finished. 
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